Puzzles in B — > h c (xc2)K Decays and QCD Factorization 
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We study the factorization-forbidden decays B —* h c K and B — * x^K in the QCD factorization 
approach. If neglecting the vertex corrections and regularizing the end-point singularities in spec- 
tator corrections properly, we get small branching ratios for both the two decay modes, which are 
roughly consistent with the experimental upper limits. This is in contrast to another factorization- 
forbidden decay B — + XcoK, for which a large decay rate is obtained in the same approach. 

PACS numbers: 



1. Introduction. The branching ratios(BRs) of B — > 
h c (xc0,Xc2)K decays are expected to be small since 
they are naively factorization- forbidden |]J. But non- 
factorizable contributions to these modes could be large 
and enhance their BRs. This is the case for B — ► XcoK, 
which was observed by Belle 2] with surprisingly large 
BR of (6.0ti;g ± 1.1) x 1CT 4 and was soon confirmed by 
BaBar Q with BR of (2.7 ± 0.7) x 1CT 4 , being compara- 
ble to those of factorization-allowed but color-suppressed 
decays, such as B^ J/ip(n c , Xci---)K. 

One may then expect similarly large BRs for B — > 
h c (xc2)K decays. Accordingly, a number of authors 
have suggested that h c could be easily observed in B- 
factories However, recent measurements set the up- 
per limit BRs of B -> XciK @ and B h c K @ at 90% 
C.L.: 

Br(B+ -> Xc2K + ) < 3.0 x 10~ 5 , 
Br(B° Xc2K°) < 4.1 x 10~ 5 , 
Bi{B+ -> Xc2K+) < 3.8 x 10~ 5 , (1) 

which are about an order of magnitude smaller than that 
of B — > XcoK. The puzzling hierarchies challenge our 
understanding of the properties of the non-factorizablc 
mechanism for these decay modes. 

The non-factorizable contributions to these modes 
have been studied in different approaches. The final- 
state soft re-scattering effects were suggested and large 
and comparable rates for all these three modes were pre- 
dicted , which are evidently inconsistent with new ex- 
perimental data. On the other hand, with the Light-Cone 
Sum Rules 8} the non-factorizable contributions were 
found to be too small to account for the large B — > XcoK 
decay rate. Recently, B — > XcoK^ were also studied in 
the PQCD approach @, where the kx factorization was 
adopted. They predicted a large rate of B — > XcoK by 
including the spectator contributions only, but no pre- 
dictions for rates of B — ► h c (Xc2)K were given. 

Within the framework of QCD factorization -B— ► 
XcO,2K decays were studied earlier |ll|,ll2jj an d we found 
that there exist infrared divergences in the vertex cor- 
rections and end-point singularities in the leading twist 



spectator corrections, which implies that soft contribu- 
tions may be large in these decays. These contributions 
come from the soft gluon exchange between the K or B 
meson and the cc pair which is emitted as a color-octet 
at the short-distance weak interaction vertex. This may 
imply a connection of the infrared behavior in the ex- 
clusive decays with the color-octet contributions in the 
inclusive B decays to charmnium , which are found to be 
the dominant mechanism [l^. The results of Ref. 01 
suggest that a sizable fraction of the large color-octet 
partial rate of inclusive B decay into charmonium does 
in fact end up in some two-body decay modes. This is 
even true for the factorization-allowed decays, such as 
B^J/tpK, for which the infrared safe leading- twist con- 
tributions only result in a BR about 1 x 10~ 4 , which 
is an order of magnitude smaller than the experimental 
data ^4|. However, after considering the soft contribu- 
tions arising from higher-twist spectator interactions, the 
prediction can account for the data quite well ^ij. The 
prominent effects of soft spectator interactions are also 
used in Ref. to explain the large BR of B -> V(3770) K 
decay. 

Similarly, these soft contributions to B — > XcoK were 
also estimated by the authors of Ref. and by us 
in QCD factorization. Both they and we found that there 
were linear singularities in chirally enhanced twist-3 spec- 
tator interactions, which are numerically large and dom- 
inant in B — * XcoK decay, but the results of theirs and 
ours were quite different because different treatments of 
the end-point singularities in the spectator interactions. 
They also got a small BR for B — > Xc2 K by introducing 
undetermined imaginary parts for soft spectator interac- 
tions. 

Since the s(s) quark emitted from the weak interaction 
vertex moves fast in the B meson rest frame, we may ex- 
pect that the soft gluon exchange is dominated by that 
between the cc pair and the spectator quark that goes 
into the kaon. For B — > XcoK d ecay , this fact has been 
confirmed in QCD factorization |16L fl7| and PQCD fac- 
torization approaches. In particular, in |17| we find 
that the vertex corrections could be small, and the spec- 
tator corrections are large and dominant. In this letter, 
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we will evaluate the BRs of B ^ h c (xc2)K by using the 
same schemes given in [I?) . We will show that if the con- 
tributions from vertex corrections are neglected, we can 
fit the experimental data of B —>h c (xc2,a)K quite well. 

2. Amplitudes in QCD factorization. We treat the 
charmonium as a non-relativistic cc bound state. Let p 
be the total momentum of the charmonium and 2q be 
the relative momentum between c and c quarks, then 
v 2 ~ Aq 2 /p 2 ~ 0.25 can be treated as a small expansion 
parameter. For the P-wave charmonium, because the 
wave function at the origin lZi(0) = 0, the amplitudes 
need to be expanded to the first order in the relative 
momentum q (see, e.g., 0]), 

M(B ^ 2S+1 Pj(cc))=J2 (IL Z ;SS Z \JJ Z ) f^q a 

x^q°)r i MTr[O a (0)PssXpfi) + mPssSPfi)l ( 2 ) 

where 0{q) represents the rest of the decay matrix ele- 
ments and are expected to be further factorized as prod- 
uct of B — * K form factors and hard kernel or as the 
convolution of a hard kernel with light-cone wave func- 
tions of B meson and K meson within QCD factorization 
approach. The spin projection operators Pss z (jPi <?) are 
constructed in terms of quark and anti-quark spinors as 
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In Eq. @ we take charmonium mass M ~ 2m c . Here 
m c is the charm quark mass. 

The integral in Eq. is proportional to the derivative 
of the P-wave wave function at the origin by 

^«m(?) = i £ * a (L*)yJ<(0), (6) 

where e a (L z ) is the polarization vector of an angular 
momentum-1 system and the value of 7^(0) for char- 
monia can be found in, e.g., Ref. fl9| . 



We use the following polarization relations respectively 
for h c (J = 1) and Xci{J = 2): 



Y / e* a (L z )(lL z ;00\U z ) = e* a (J z ), 

LzSz 

^ e * Q (L,)e*^ 2 )(lL 2 ;lS,|2J z ) = e*^(J z ). (7) 

LzSz 

where s a (J z ) is the polarization vector for h c , and the 
polarization tensor e al ^(J z ), which is symmetric under 
the exchange a «-> (3, is the one appropriate for \ci ■ 

The effective Hamiltonian relevant for B — > h c (Xc2)K 
is written as 201 



n cS =^(v cb v; s (c 1 o 1 + c 2 o 2 )- v tb v t 



(8) 



where Gf is the Fermi constant, and C,-'s are the Wil- 
son coefficients. The relevant operators CVs are given in 
Ref. [13. 

According to ^(j an non-factorizable corrections are 
due to the diagrams in Fig^ These corrections, with 
operators Oi inserted, contribute to the amplitude 0{q) 
in , where the external lines of charm and anti-charm 
quarks have been truncated. Taking non-factorizable cor- 
rections in Fig^into account, the decay amplitudes for 
B -^h c (xc0.2)K in QCD factorization can be written as 
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(9) 



where in the parentheses, "— " is corresponding to h c and 
"+" to Xc2, respectively. Here and afterwards, we will use 
the subindex j to denote different charmonium states, 
and j = 1, 2, represent h c and Xc2 respectively. Thus the 
coefficients A/s are given by 

4 -^-^(^(rt +/*,).<«>) 

Here F\ is one of the B^K form factors and the relation 
Fo(Mf)/Fi(Mf) ~l-M?/m% has been used to sim- 
plify the expressions of the amplitudes. The function flj 
is calculated from the four vertex diagrams (a, b, c, d) 
and fllj is calculated from the two spectator diagrams 
(e, f) in FigQ] The function fllj receives contributions 
from both twist-2 and twist-3 light-cone distribution am- 
plitudes of the K meson, and we will simply symbolize 
them as fll 2 and f Uj, respectively. 
The factors rij in (|1U|) are given by 



m 

ii 2 



(1 - z 2 )m B 



(11) 



where z n - = M 2 /m% 



For B — > h c {Xc2)K decays, just like B — > XcqK, the 
soft gluons couple to the charm quark pair through color 
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FIG. 1: Feynman diagrams for vertex and spectator correc- 
tions to B —> h c (xc2)K . 



dipole interactions, which are proportional to the relative 
momentum q and give the leading order contributions 
both in l/m& and NR expansion (see Eq. (J2J). As a 
result, there are generally infrared divergences in flj, 
while fllj ' 3 suffer from logarithmic and linear end-point 
singularities, respectively. The infrared divergences in 
flj can be regularized by a gluon mass m g or by the 
binding energies bj = Mj — 2m c 16] of charmonia, 



fh = 



iz(l — z+hiz) f m 



{1-zf 



In ( — | ) + finite terms, 



4z((l-z)(l+3z)+4zlnz) 



In — finite terms, (12) 

1V1 rj 



which are consistent with Ref . [l2l Il7j and Ref. ^| , re- 
spectively. Here and afterwards we will omit the sub- 
script of Zj unless it is necessary. 

To derive the functions fllj , we use the light-cone pro- 
jector of K meson up to twist-3, 



4 



\i>' i$4>K{y) 



9 r K (v) 



dk 



2// 



G 



')} fl .(13) 



where fc 2 (i) is the momentum of the anti-quark (quark) in 
K meson, and the derivative acts on the hard-scattering 
amplitudes in the momentum space only. The chirally 
enhanced mass scale hk = mK 2 / (m s (/x) +md(n)) is com- 
parable to rrib, which ensures that the twist-3 spectator 
interactions are numerically large, though they are sup- 



pressed by 1/mf,. 1 The twist- 2 light-cone distribution 
amplitude (LCDA) (j>K{y) and the twist-3 ones 4? K {y) 
and 4>j((y) are symmetric under exchange y<-> (1 — y) in 
the limit of SU(3) symmetry. In practice, we choose the 
asymptotic form (f> K (y) = <f>% : (y) = 6y(l-y) and (f> p K (y) = l 
for these LCDAs. 

The projectors in Eq. I|13|l have been used by us in 
Ref. |l7| and there exist logarithmic and linear end-point 
singularities in twist-2 and twist-3 spectator interactions, 
respectively. Similar singularities are also found by the 
authors of Ref. Hfl- But they use a different projector, 
which can be derived from Eq. (|13l) by adopting an inte- 
gration by parts on y and dropping the boundary terms. 
As having mentioned in Ref. [13, the two projectors can 
be consistent with each other only when the singularities 
in all the regions with small-virtualities are carefully reg- 
ularized. That is, to introduce a relative off-shellness A 
to regularize every factor "y" in the denominator. For 
linear singularities, we have 



y n dy 



1 



(y + A) 



n+2 



(n + l)A 



-1 + 0(A), n = 0,l,2... (14) 



The difference between Eq. II14|I in our scheme for small- 
virtuality and that in Ref. [16( is quite evident. If simply 
parameterizing the linear singularities on the left hand 
side of Eq. ifHj) as f dy /y 2 = tub/ Ah following Ref. [T^ |. 
one would get the same results for all n > 0. We will see 
that the non-trivial factors (n + 1) _1 in Eq. (|14|> play an 
important role in evaluating the BRs of B — ► h c (x C 2)K 
decays. 

Under this scheme , as we expected, the projectors in 
Eq. i|13|) and in Ref. [l|j give the same forms of fUj' '■ 



/2Z? =or ^[3* + 0(A)], 

///f = a 2 -^-[6zlnA + 15.?+ 0(A)], 



3 a^TK ™>b r — 3z 



(15) 



m = ^- z ~£[— -(l + ,)lnA--(l-,) + 0(A)], 



m 



3 _02Tk_ m B r-z 
~ 1 - z As T " " '2 



(l+2z) lnA--(3 + 2z) + 0(A)](16) 



where tk = 2/xjf/mt, is of 0(1) and the factors aj are 
defined as 



^ 2 !k!b 



(17) 



J iV c (i- 2i )2 m | Fl (M?)' 

In Eq. (USJl and Eq. JTHJ), the integral with LCDA of B 
meson is conventionally parameterized as dt; = 
fa with A B « 300 MeV 



1 In fact, the twist-3 contributions to these decay modes are not 
power suppressed because of the linear singularities contained in 
them. 
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TABLE I: Theoretical predictions for functions fllj and Br,. 
The values given at /i = 1.45 Gev (u = 4.4 Gev) are shown 
without (with) parentheses. 





10 4 x Brj 






j 


m g = 0.5 Gev 


m g = 0.2 Gev 






B^h c K 


1.6(0.8) 


2.7(1.3) 


3.1 


-12.9(-18.3) 




0.8(0.3) 


1.4(0.5) 


2.9 


-5.7(-8.7) 


B^XcoK 


3.0(1.9) 


4.2(2.4) 


7.8 


-35.0(-53.6) 



It is worth emphasizing that this scheme is physical 
since the off-shellness of quarks and gluons are naturally 
serve as infrared cutoffs when Following Ref. |17j . 

we use the binding energy bj to determine the relative 
off-shellness A, : 



M-i 



(18) 



Before explicitly evaluating these amplitudes, we 
should emphasize the importance of the twist-3 specta- 
tor contributions again. Normalizing all the amplitudes 
by the finite part of flj , we can see that both flj 's and 
fUj's are of the order In A ~ ln(mB/A), but fUj's are 
of order -J^- ~ ^fi- because of the chiral enhancement 

Am b A 

and the linear singularities contained in them. These 
unusual power counting rules give support to the naive 
expectation that the soft spectator interactions would be 
dominant in these three decay modes. At the qualitative 
level, this statement has been validated in Ref. 0, 0] 
for B— >XcoK- 

3. Numerical results and discussions. For numerical 
analysis, we use the following input parameters with two 
values for the QCD renormalization scale \i = 1.45 Gev 
and £i = 4.4 Gev: 

Mi = 3.524 GeV, M 2 = 3.556 GeV, m c = 1.5 GeV, 
m B = 5.28 GeV, A B = 300 MeV, n[ (0) = V0.075 GeV 5/2 , 
./s = 216 MeV 21], f K = 160 MeV, 



Fi(Ml) = 0.80, Fi(Mf) = 0.82 [23, 
r K (n) = 0.85(1.3), a s (/j,) = 0.34(0.22). 



(19) 



In Eq. l(T§j l the ^-dependent quantities at fi = 1.45 Gev 
(/i = 4.4 Gev) are shown without (with) parentheses. The 
Wilson coefficients Cj's are evaluated at leading order 
by renormalization group approach |20j , since the ampli- 
tudes in Eq. (|10fl are only of the leading order in a s . 

We use LoopTools j23| for numerical calculations of 
flj with a gluon mass m g , and the obtained BRs Br^ 
are listed in Tab. I with m g varying from 200 Mev to 500 
Mev. The values of fllj are also shown in Tab. I. We see 
the ^-dependence of f JJ| arising from Tk{^) is largely 
cancelled by a s (ji) [T(j. In Tab. I, for comparison, we also 
list the results of B^\ c qK decay, which can be found in 



TABLE II: Predictions of Brj with only spectator contribu- 
tions fllj. 



10 4 x Brj 




= 4.4 Gev 


/' 


= 1.45 Gev 


B -> h c K 




0.27 




0.27 


B X c2K 




0.03 




0.02 


B -» XcoK 




1.05 




1.17 



Ref. 0| 2 - Here and afterward, we will use subscript "h c , 
Xc2, Xco" to symbolize the quantities for h c , Xc2 and XcO 
respectively. 

We predict a relative small BR for B — > Xc2K with 



Ry 



Br Xe / Br ?< 



3 — 6 in Tab. I, which gives a 



signal of hierarchy without fine tuning of parameters. 
However, the BR of B — > h c K is large and the ratio 
Ri = Br Xc0 /Br/ lc ~ 2, which is not large enough to ac- 
count for experimental data. 

Further, we regularize the infrared divergences in 
flj(j = h c ,Xc2,Xco) using the binding energy. The re- 
sults for fI Xe3 and fI Xc0 are available in Ref. [l^. We 
find both fI Xa2 and fI Xc0 in the binding energy scheme 
are smaller than those in the gluon mass scheme (with 
the gluon mass values chosen above). Since the large 
Br Xe0 are mainly due to the large twist-3 spectator con- 
tribution fH Xc0 0, it is not very sensitive to the value 
of fl, but Br Xc2 is. As a result, we get R Xc2 > 20 no 
matter which fj, is chosen. Numerical calculations indi- 
cate that the value of /J/ lc in the binding energy scheme 
is also smaller than that in the gluon mass scheme. 

The largeness of R Xc2 is mainly due to the difference 
of fll^ and fZT? , and can be traced to the infrared 
behavior of the spectator interactions. This is the reason 
why our results are different from those of Ref. |16| . in 
which the coefficient of the linear pole in fII Xo2 is two 
times larger than ours in (|16|l . As a result, the ratio R Xc2 
can be large only by adding adjusted imaginary parts to 
spectator functions fllj there. 

As we have mentioned, since the s(s) quark going out 
from the weak interaction vertex moves fast in the B me- 
son rest frame, the soft gluon emitted from the cc quark 
pair is favored to be absorbed by the spectator quark. 
This picture has been confirmed by the power counting 
rules deduced above and by our analysis for B—>XcoK de- 
cay ^3 • On the other hand, the vertex functions flj 's in 
general are numerically small but sensitive to regulariza- 
tion schemes. So it might be reasonable to evaluate am- 
plitudes with only spectator contributions fll? and fllj, 
and treat vertex corrections as uncertainties, of which 
the contributions to Br^ are expected to be smaller than 
1 x 10~ 4 , as in the case of factorizable decay B — > J/tjj K. 
And it could be even smaller for nonfactorizable decays. 



2 3 rrm 

the functions fII Xc0 are extracted from Eq. 1101 with 

, u = — (1 — z)rriB /v / 3z, so their definitions differ from those in 
Ref. Il7j just by a minus sign. 
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With the approximation of neglecting contributions from 
vertex corrections, the results are listed in Tab. II. 

In this approximation, the soft spectator interactions 
dominate, and there is no surprise that we get a larger 
Br Xc0 and a much smaller Br Xe2 simultaneously. The 
hierarchy is more evident and R Xa2 > 30. Furthermore, 
the rate of h c is smaller than that in Tab. I, and J?/ to ~4, 
which is roughly consistent with the experimental data if 
we take into account the new PDG average value Br Xc0 = 

(1.6+2:5) x iQ- 4 - 

Since soft scattering can introduce a strong-interaction 
phase, the functions fllj could be complex. In our 
scheme, the phases may emerge with logarithmic singu- 
larities in fllj if A is negative. To account for this effect, 
here every logarithm in l|15f) and H16|) can be multiplied 
by a universal factor (1 + pe 10 ) [l(j with < p < 1.5. 
The phase 9 is expected to be small, say, varying from 
to 7r/2. We find that the corrections to those BR.^s 
do not change our conclusions drastically. For example, 
choosing p — 1.3 with 6 varying from to 7r/2, we get 
Br XcO ~(0.4-4.0)xl(T 5 , i? Xc2 ~100 and R ha ~A, which 
are roughly consistent with experimental data, respec- 
tively. 

In summary, since the non-factorizable contributions 
arising from soft spectator interactions could be large for 
B to charmonium exclusive decays, the smallness of ex- 
perimental upper limits of BRs of B^h c {x C 2)K are sur- 
prising. We study these two modes within the framework 
of QCD factorization. If neglecting the vertex corrections 
and regularizing the end-point singularities in spectator 
corrections properly, we get small BRs of 0.27 x 10~ 4 



and 0.03 x 10 -4 for B — > h c K and B — > Xc2K respec- 
tively, while the predicted BR for B—*Xct)K is large, as 
shown in Table II. These are roughly consistent with ex- 
periments. On the other hand, at present we do not have 
a rigorous treatment for the vertex corrections because 
of the appearance of infrared divergences. Nevertheless, 
if we use the gluon mass (or binding energy) to regular- 
ize the infrared divergences (with reasonable values for 
the gluon mass or binding energy), we find their effects 
are possibly small, and may therefore be viewed as the- 
oretical uncertainties. If we take these uncertainties into 
account, the predictions for BRs of B — > h c (Xc2)K be- 
come Br hc = (3±10)xl0~ 5 and Br Xc2 = (0±10)xl0~ 5 . 
Here the large errors in predictions are mainly attributed 
to uncertain but sub-important (even power-suppressed) 
contributions arising from those vertex corrections. 

Note. Very recently, the decay mode B^h c K has also 
been studied in PQCD factorization approach |2f|. As in 
the case of B— >XcoK .9], the authors of [2{| neglected the 
vertex corrections, and get a small rate Br(£>^ h c K) = 
4.5 x 10~ 5 , which is roughly consistent with ours. 
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